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ABSTRACT 



Aims. We consider isolated compact remnants (ICoRs), i.e. neutrons stars and black holes that do not reside in binary systems and 
therefore cannot be detected as X-ray binaries. ICoRs may represent ~ 5 percent of the stellar mass budget of the Galaxy, but they 
are very hard to detect. Here we explore the possibility of using microlensing to identify ICoRs. 

Methods. In a previous paper we described a simulation of neutron star evolution in phase space in the Galaxy, taking into account 
the distribution of the progenitors and the kick at formation. Here we first reconsider the evolution and distribution of neutron stars 
and black holes adding a bulge component. From the new distributions we calculate the microlensing optical depth, event rate and 
distribution of event time scales, comparing and contrasting the case of ICoRs and "normal stars". 

Results. We find that the contribution of remnants to optical depth is slightly lower than without kinematics, owing to the evaporation 
from the Galaxy. On the other hand, the relative contribution to the rate of events is a factor ~ 5 higher. In all, ~ 6 - 7 percent of 
the events are likely related to ICoRs. In particular, ~ 30-40 percent of the events with duration > 100 days are possibly related to 
black holes. 

Conclusions. It seems therefore that microlensing observations are a suitable tool to probe the population of Galactic ICoRs. 
Key words, stars: kinematics - stars: neutron - stars: statistics 



1. Introduction 

By isolated compact remnants we mean neutron stars and black 
holes (NSs and BHs hereafter) that are generated from the core- 
collapse of isolated massive stars. Hence, these objects cannot be 
recycled by accretion of matter stripped from a companion star. 
For NSs we consider in particular those objects that exhausted 
their possible reserves of energy, like rotation and internal heat, 
and are not detectable as active sources, as in the case of radio 
pulsars or X-ray pu lsars, like the Magnificent Seven (e.g. Haberl 
I20071 lTurollall2009l and references therein) . 

Extensive searches for ICoRs as sour ces powered by ac- 
cretio n from the inte r stellar m edium (e.g. IPopov & Prokhorov 
1 19981 iTreves et alj l2000b lAgol & Kamionkowskil |2002|) 
were performed, for examp le, in the ROSAT catalog (e.g. 
iNeuhauser & Trumpet! 1 19990 . However, no positive confirma- 
tion of an accr e ting IC oR has been reported to date. 

lOfek et aTl d2010h suggested that the recently discovered 
long duration radio transients (e.g. Bow er et al.l2007l) may be as- 
sociated with isolated old neutron stars (ONS). Hence, a possible 
direct detection of ICoRs may already be available, but is n ot yet 
fully demonstrated. On the other hand lBennett et all (12002) sug- 
gested that ICoRs, in particular black holes, may be responsible 
for some of the long-duration microlensing e vents observed. 

Microlensing surveys (see M oniejl2010l for a recent review) 
are used to study the distribution of dark as well as luminous 
matter along a given line of sight (l.o.s. hereafter), and there- 
fore this may be a suitable method to probe the phase space dis- 
tribution of ICoRs. The first theoretical studies of the subject 
(iGriestl 1 1 99ll iPacz vnskil 1 1 99 lh predicted an optical depth t ~ 
0.5-0.8 x 10~ 6 for the luminous stellar component of the Milky 
Way (see Sect. [2] for a definition of t). Several observational 
campaigns have been set up to detect microlensing events in 



the di rection of the Galactic bulge , e.g. MA CHO fAlcock et all 
fl997l) . OGLE (lUdalskietalJ 19941) and MOA (ISumi et alJ2003l) . 

The first measurements of the optical depth returned values 
that a re significantly highe r than those predicted, e .gr = 3.3 x 
1Q- 6 dUdalski et alJ[l994h and r = 3.9 x lO" 6 (lAlcocketal 



1 19971) . To explain this discrepancy. iKiraga & Paczvnskil 1994 
suggested that the Galactic bulge could give a substantial con- 
tribution to the optical depth and thus they included this contri- 
bution in their calculations. Theoretical estimates with these new 
models agree well with the m ore recent measurem ents of the op- 
tical depth, t = 2.17 x 10~ 6 (Pop owski et al.l2005l the MACHO 
group), t = 2.55 x 10' 6 (ISumi et al.l2006l th7o GLE-2 group) 
and t = 1.62 x 10~ 6 (Hamad ache etaD 120061 the EROS-2 
group). 

The contribution o f comp a ct remnants to micr o lensin g 
has been estimated bv iGouldl (120001) : IWood & Maol J2005); 
ICalchi Novati et al l (120081) . These authors found that ~ 2 - 
3 percent of the events are related to ICoRs. However, they 
did not take into account that ICoRs have a different phase 
space distribution with respect to normal st ars, because of the 
high kick velocities re ceived at birth, (e.g. iHobbs et alj [2005; 
Gualandris et al. 2005). This has two opposite effects. First, the 
velocities of ICoRs at birth can be higher than the escape veloc- 
ity from the Galaxy. Hence, a non negligible fraction of remnants 
may have escaped from the Galaxy and thus cannot contribute to 
microlensing toward the bulge. Second, the higher velocities of 
ICoRs imply an event rate higher than that expect ed from a sim - 
ilar population of lenses moving at lower speeds (rGriestll 19911) . 
Thus, in this paper we investigate the possibility to probe iso- 
lated compact remnants by means of microlensing observations. 
We base our work on the results of Monte Carlo simulations of 
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NS orbits (IOfekll2009l [Sartore et alJl2010l hereafter Of09 and 
Paper I, respectively). 

The paper is organized as follows. In Sect. [2] we describe 
our models for the distribution of bulge and disk stars as well as 
ICoRs after an introduction of the basic microlensing quantities 
and expressions. We present our results in Sect. [3] In particular 
we compare the optical depth and event rate due to ICoRs with 
that of normal stars. We also study the distribution of event time 
scales in both cases. Finally we discuss our results in Sect. [4] 

2. Model 

We now introduce some basic expressions related to the mi- 
crolensing phenomenon. The probability that a background 
source is lensed at any given time is called microlensing optical 
depth and depends on the distribution of lensing matter along the 
l.o.s. (e.g. Kira ea & Paczvnskilll994lJetzer et al.ll2002l) 



t(D s ) 
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(1) 



where pi is the mass density of the lenses, D s is the source- 
observer distance, and D\ is the lens-observer distance. If both 
lenses and sources are distributed along the l.o.s., Eq.[T]becomes 
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where p s is the mass density of sources and / = 
J () p s (D s )D s dD s is a normalization factor. Throughout this 
paper we assume D max = 12 kpc to include the contribution of 
the whole bulge (see S ect. 12. U . The optical depth can also be 
defined as dGriesjl 19911) 



T = 



n (t E ) r 



(3) 



where F is the rate of lenses entering the tube and (fg) is the av- 
erage time scale of the observed events. The duration of a single 
event depends on the mass of the lens and the geometry of the 
system 
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where v± is the relative velocity between source and lens in a 
plan e perpendicu lar to the l.o.s. and Re is the Einstein radius 
(e.g. lGriesdfT99l 
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1.38 x 10" 8 
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where M is the mass of the lens, x = D//D s and m = M / M 
is the mass of the lens in solar units. 

The differential rate of events is (e.g. iJetzer et all [2002, 
ICalchi Novati et al.ll2008l) 



ni (Di) d 3 D, n s (D s ) D 2 dD s - 

dr = f(\±)d z v ± . 

dt / 



(6) 



where «/ and n s are the number density of the lenses and the 
sources along the l.o.s., while /(Vj_) is the distribution of the 
source-lens relative velocities. 



2.1. Distribution of bulge and disk stars 

We model the Galactic bulge as in ICalchi Novati et a D (|2008), 
i.e. as a triaxial bulge with an exponential density profile and a 
major axis that forms an angle = 23.8° with the Sun-Galactic 
center axis (e.g lStanek et al.l ll997). The total mass of the bulge 
is ~ 2 x 1O 1O M . For the disk, we adopt a thin + thick disk 
model, both described by an exponential density profile 



Pd,(R,z) = 
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where the masses of the thin and thick disks are Mo, ~ 2.5 x 
10 10 M o andM D2 ~ 0.5 x 10 10 M Q , the scale-lengths are L x ~ 
2.6 kpc and L2 ~ 3.6 kpc and scale-heights are H\ ~ 0.3 and 
H2 ~ 0.9 kpc (e.g. ?). The mass of the thin disk accounts also 
fo r the interstellar med ium, Mjsm ~ 0.5 x 1O 1O M . According 
to iFreudenreichl (1 1 9981) . the stellar disk is holed at its center, the 
hole being likely produced by orbital resonances in the potential 
of the barred bulge (e.g. lConto poulos et al. 1989). For the scale- 
length of the hole of both the thin an d thick disks we assume 
L h ~ 1 .3 kpc dPicaud & Robinll2004l) . 

The motion of bulge and disk stars has both bulk and random 
components. The Galactic bulge does not rotate as a whole like 
a rigid body but out of a certain radius the bulk velocity flattens 
dRich et al. 2007). Thus we assume that the rotation velocity v r J t 
of the bulge grows linearly to 50 km s~ out of a radius of 1 kpc 
from the Galactic center. Out of this radius, we assume a flat ro- 



tation curve, with 



„«0 



50 km s . For disk stars, we compute 



the bulk motion self-consistently from the potential generated 
by disk and bulge stars, see Sect. 12.21 We assume for simplicity 
that the random motions of bulge and disk stars are isotropic, 
with dispersions cr^ = 100 and crj rf) = 25kms _I which rea- 
sonably agrees with the values inferred from observations (e.g. 
ICalchi Novati et al.1 l2008 and references therein). 



2.2. Isolated compact remnants 

In Paper I we studied the dynamics of a population of disk NSs 
born at a constant rate during the Milky Way lifetime, assum- 
ing a total of 10 9 NSs, consi stently with chemi cal abundances 
observed in the Galaxy (e.g. lArnett et al.l ll989). As Of09 has 
pointed out, the total number of disk-born NSs inferred from the 
present-day supernova ra te and from the star-for mation history 
of the disk is < 4 x 10 8 dKeane & Kramerll2008l and references 
therein). To explain this discrepancy, Of09 suggested that the re- 
maining NSs have been generated in the bulge. 

Because we are dealing with the microlensing rate of com- 
pact remnants toward the Galactic bulge, it is straightforward 
to think that the major contribution comes from bulge-born ob- 
jects. Thus we re-run our simulations, taking into account the 
contribution of the bulge. First, we estimate the t otal number o f 
ICoRs born in the MW, following the approach of lGouldH gOOO). 
We ad opt the initial mass function (IMF) proposed by Kroupa 
(12001b . i.e. a triple power law model 



dN 
dm 



a = 0.3 , 
a = 1.3, 
a = 2.3, 



0.03 < m < 0.08 
0.08 < m < 0.50 
0.50 < m < 100 . 



(8) 



From the IMF we estimate the number and mass fractions 
of each stellar population - brown dwarfs, main sequence stars, 



N. Sartore, A. Treves: Probing isolated compact remnants with microlensing 



3 



white dwarfs (BDs, MSs and WDs hereafter), NSs and BHs. We 
assume that all stars with m > 1 have evolved through the rem- 
nant phase, i.e. stars with 1 < m < 8 are now white dwarfs 
(mwD = 0.6), while stars with masses 8 < m < 40 and 
40 < m < 100 are treated as neutron stars (m^s = 1.4) and 
black holes {mbh = 10), respectively. Results are reported in 
Table [1] We can now obtain the number of NSs and BHs gen- 
erated in the Galaxy: ~ 3.3 x 10 8 , N% } H ~ 3.2 x 10 7 , 

N { £, ~ 4.1 x 10 8 and ~ 4.0 x 10 7 , where the superscripts 
(b) and (d) refer to the bulge and disk populations, respectively. 
We note that the assumption that all disk stars with 1 < m < 8 
are now white dwarfs gives only a rough estimate of the number 
of disk WDs. On the other hand the assumption is well justified 
for NSs and BHs, because the typical lifetime of massive stars is 
much shorter than the age of the Galaxy. 



BD 


MS 


WD 


NS 


BH 


number fraction 0.272 
mass fraction 0.059 


0.653 
0.744 


0.065 
0.157 


0.004 
0.023 


0.0004 
0.016 



Table 1. Number and mass fractions of the different stellar pop- 
ulations. 



We use the PS YCO code (Paper I) to follow the orbits of 2 x 
10 s bulge-born and 2 x 10 s disk-born synthetic ICoRs. We have 
revised the Galactic potential assumed in P aper I, considering a 
super position of Miyamoto and Nagai disks (Miy amoto & Nagail 
1 19751) for both the bulge and disk. The associated density pro- 
files are a close approximation of the density profiles assumed 
for the bulge and disk of the Galaxy. As in Paper I, we adjusted 
the parameters of the dark matter potential to obtain a circular 
velocity of ~ 220 km s _1 at the solar circle (R Q = 8.5 kpc). The 
resulting escape velocity at the same radius is ~ 450kms~ 1 , i.e. 
~ 17 percent lower than in Paper I. Also, we make the simpli- 
fying assumption that the gravitational potential of the bulge is 
axisymmetric. 

The results obtained are then normalized to the correspond- 
ing number of objects expected from each population and re- 
ported above. The initial conditions for each object are randomly 
assigned using a Monte Carlo procedure. The integration time 
for each object is equal to its age, i.e. 10 Gyr for all bulge ICoRs, 
while the age of disk ICoRs is uniformly distributed between 
and 10 Gyr. In both cases the radial birth probability is assumed 
to be proportional to the density of stars, and we assume also 
that the bulge and disk can produce remnants up to 3 and 15 
kpc, respectively. For the bulge, the vertical birth probability is 
again assumed proportional to the density of normal stars up to 3 
kpc, while f or the disk we ado pt the same distribution as the one 
proposed bv lPaczvnskil(ll990i) . which follows the distribution of 
young stars in the disk of the MW. The azimuthal coordinates 
are uniformly distributed in both cases. 

The initial velocity of each synthetic remnant is the vector 
sum of its birth velocity plus the orbital velocity of the progeni- 
tor. The orbital velocity of the progenitor is calculated from the 
gravitational potential. The birth velocities of NSs are distributed 
according t o a Maxwell ian distribution with dispersion cr 
265 km s dHobbs et all 120051) . The velocity distribution of 
BHs is poorly constrained, and it was assum ed that these objects 
have a dispersion of ~ 40km s _1 (e.g. IWhite & van"p"ar adiis 
Il996l) . However, the discovery of BH X-ray bi naries with high 
spatial velocities (e.g. lMirabel et al.l200lll2002h . points to a sim- 
ilar birth velocity distribution for NSs and BHs. Hence, in this 
work we assume that BHs have the same velocity distribution of 



NSs. Finally, we also add a random component to the velocity of 
the progenitors, see Sect. 12. II For the bulge, this dispersive com- 
ponent is dominant with respect to the orbital velocity around 
the Galactic center. 

3. Results 

3. 1 . Neutron stars 

First, we compare the statistical properties of NSs in the new 
model with those reported in Paper I and Of09. We find that for 
bulge NSs, the evaporation is highly inefficient, that is, almost all 
NSs are in bound orbits, f bound ~ 0.98, because of the high grav- 
itational force toward the inner part of the Galaxy. The fraction 
of disk-born NSs in bound orbits is slightly lower than that found 
in Paper I for the same velocity distribution, fbound ~ 0.63 be- 
cause, as we already pointed out, the escape velocity in the new 
model is also lower. In all, ~ 20 percent of NSs evaporated from 
the Galaxy. The density of bulge-born NSs becomes dominant in 
the inner part of the Galaxy, R < 3 kpc. In particular the density 
of ONS in the bulge can be as high as ~ 0.1 pc~ 3 , see Fig.Q](left 
panel). 

Our results are different from those obtained by Of09, for 
which fbound varies from 0.41 to 0.52 for bulge NSs and from 
0.13 to 0.16 for disk NSs, respectively, depending on the initial 
velocity distribution. The main reason is the difference in the 
gravitational potential adopted for the M W. In our model the to- 
tal mass of the Galaxy is ~ 10 12 M (e.g. lXue etalll2008h . i.e. a 
factor ~ 7 larger than that adopted by Of09. This obviously im- 
plies that in our model a larger fraction of ONS is bound to the 
Galaxy. Furthermore, we adopt different initial conditions, e.g. 
the distribution of progenitors and birth velocities of NSs and a 
different (lower) normalization. 

The local density of NSs is «o ~ 3.3 x 10 pc~ 3 , of which 
the fraction of bulge-born ONS is ~ 0.23. The density is there- 
fore a factor ~ 3 lower than what we found in Paper I for the 
same velocity distribution and a factor ~ 10 lower than that re- 
ported by Of09. We note however that in Paper I we obtained 
a value of the local density similar to that of Of09 when using 
the same velocity distribution. Also, the fraction of bulge-born 
objects is similar to that of Of09. This means that local density 
is affected mainly by the normalization, i.e. the total number of 
NSs in the Galaxy. 

We find that the density in the Galactic center is «cc ~ 
0.12pc~ 3 , of which ~ 93 percent have been born in the bulge 
itself. These results are similar to those of Of09, who obtained 
n GC ~ 0.2 - 0.3 pc~ 3 , with 95 - 97 percent of NSs being 
bulge-born. The larger contribution of bulge NSs to the density 
toward the Galactic center can be appreciated for example in 
Fig. [2] where the projected density within 12 kpc from the Sun 
is ~ 7.3 x 10 4 deg~ 2 against ~ 6.7 x 10 3 deg~ 2 of disk-born 
NSs. 

3.2. Black holes 

Because we assumed that NSs and BHs have the same kinematic 
properties, the density ratio between BHs and NSs should be 
~ 0.1, as expected from the ratio between the number fractions 
of the two populations (e.g. Table [TJ. Indeed we find that the 
local density of BHs is ~ 3.2 x 10~ 6 pc~ 3 , while the density of 
BHs in the Galactic center is ~ 1.1 x 10~ 2 pc~ 3 . The fraction of 
bulge-born objects is the same of NSs. The projected density of 
BHs in all directions is also ~ 10 times lower than the projected 
density of NSs (Fig. [2] lower panels). 
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Fig. 1. Density profiles of NSs (left) and BHs (right) in the Galactic plane. Solid, dotted, and dashed lines represent the total, bulge, 
and disk contributions. 





Fig. 2. Aitoff projection of the sky density of NSs (upper panels) and BHs (lower panels) within 12 kpc from the Sun, for bulge-born 
(left panels) and disk-born (right panels) objects. Contour profiles are drawn at (10, 10 2 , 10 3 , 10 4 ) deg~ 2 , respectively. Coordinates 
represent Galactic longitude (abscissae) and latitude (ordinates). 



3.3. Optical depth 



We estimate the microlensing optical depth (Eq. |2]i for differ- 
ent l.o.s. toward the Galactic bulge. In particular we consider a 
20° x 20° window, centered at (/, b) = (0°, 0°) and study the 
dependence of r on the Galactic longitude and latitude. We as- 
sume 100 percent detection efficiency in our calculations and do 
not take into account the effect of interstellar extinction nor do 
we apply any flux limit on source stars. In this regard, we con- 
sider as sources only the stars in the main sequence, because we 



expect brown and white dwarfs to be too weak to be efficiently 
monitored by present surveys. 

In Fig. |3] we show the contour plots of the optical depth 
of normal stars and ICoRs. For normal stars the optical depth 
mostly depends on the Galactic latitude. On the other hand, the 
dependence of T star from the longitude is we aker (compare this 
figure with Fig. 1 of lCalchi Novati et al.l2008l) . The optical depth 
of ICoRs instead has a strong dependence on both the longitude 
and latitude. Also, the asymmetry with respect to the Galactic 
center is more noticeable (Fig. [3] right panel). The maximum of 
the optical depth is shifted toward negative longitudes, I ~ -2°. 
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This behavior is due to the l.o.s. intercepting the far end of the 
Galactic bar when looking toward negative longitudes. In this 
condition the geometry of the source-lens system favors a larger 
Einstein radius, which implies a higher optical depth. 



l.o.s. 

(lb) 


[x ltr 6 ] 


' NS 

[x 10- 8 ] 


1 NS 

[x 10- 8 ] 


BH 

[x 10- 8 ] 


BH 

[x 10- 8 ] 


(l%-3.9°) 


0.94 


1.51 


0.31 


1.05 


0.21 


(l°.5,-2.68°) 


1.43 


1.97 


0.38 


1.36 


0.27 


(l°.5,-2.75)° 


1.40 


1.93 


0.37 


1.33 


0.26 



Table 2. Optical depth of normal stars, NSs and BHs toward 
different l.o.s. The superscripts (b) and (d) refer to the subpopu- 
lations of bulge and disk objects, respectively. 



In Table [2] we report the optical depth of normal stars, NSs 
and BHs for several specific l.o.s in the selected window. We 
find that the optical depth of normal stars is, for example, T star ~ 
0.94 x 1(T 6 toward the Baade's Window, (l,b) = (l°,-3°.9) 
while we obtain 1.43 x 1(T 6 and 1.40 x 10~ 6 toward (l,b) = 
(l°.50,-2°.68) and (l,b) = (1°.50, -2°.75). These values agree 
reasonably well with those found in the literature, see Sect. [1] 

Without dynamical evolution, the contribution of ICoRs to 
the optical depth would be proportional to the mass fraction of 
these objects. Also, it would not depend on the line of sight. 
From our assumptions on the IMF, the optical depth of NSs and 
BHs would be ~ 2.3 and ~ 1 .6 percent of the total. 

When the kinematics is taken into account, the optical depth 
differs from the value expected without dynamical evolution. We 
find that the optical depth is slightly lower than the nominal 
value obtained when the kinematics is not accounted for. Also, 
we find that the contribution of ICoRs to the optical depth is 
mostly due to objects born in the bulge, T {h) / (T (d) + t (I,) ) ~ 0.84, 
where r (i) and r (rf) are the optical depths of bulge-born and disk- 
born ICoRs, respectively. The contribution of BHs is a factor 
~ 1 .5 lower of that of NSs. 

3.4. Event rates 

We calculated the expected distribution of event time scales for 
various l.o.s. toward the bulge. Results for the Galactic center 
and the Baade's Window are reported in Table [3] as examples. 
In general, microlensing events are dominated by self-lensing of 
bulge low mass stars, these events having a typical duration of 
~ 15-20 days (see e.g. Fig . |4ji . As expected the rate of events 
drops away from the plane (see Table 01. The relative contribu- 
tion of compact remnants is increased by a factor ~ 5 with re- 
spect to the case with no kinematic effects. Indeed, we find that 
the overall contribution of NSs rises from ~ 1 to ~ 5 percent, 
while for BHs it rises from ~ 0.2 to ~ 1 percent (Fig. [5]). 

The average duration of the events associated to remnants is 
lower by a factor ~ 1.5, owing to the high velocities of these 
objects. For NSs (tg) ~ 25 days instead of ~ 36 days, while for 
black holes we find (f^) ~ 67 instead of ~ 95 days. Intriguingly, 
the relative contribution of BHs increases with the time scale. 
For events with duration longer than 100 days, this contribution 
is ~ 40 percent (see Fig. |5j, while NSs account for ~ 10 percent 
of the events. 



remnants (NSs and BHs) to the microlensing optical depth and 
event rate toward the Galactic bulge. We first ran new simula- 
tions of remnant orbits in a refined model of the Galactic grav- 
itational potential, which accounts for the observed distribution 
of stars obtained from IR surveys. In our model, the Galactic 
bulge accounts for ~ 45 percent of the compact remnants born 
in the MW. Our results show that the initial conditions play a 
fundamental role in the final fate of an ICoR. At variance with 
disk-born objects, almost all of those born in the bulge are likely 
retained by the Galaxy, and thus they give the major contribution 
to the microlensing rate toward the bulge itself. 

We found that the microlensing optical depth of ICoRs is 
lower than that obtained without dynamical evolution (i.e. no 
kick at birth). On the other hand, the net effect of the high ve- 
locities on microlensing is an increase of the event rate and a de- 
crease in the typical time scale of the events. Indeed, we found 
that the contribution to the event rate is higher than without kine- 
matics effects. The average time scale of events associated to 
ICoRs is also lower than that expected from a population of 
remnants without kicks. Interestingly, for durations > 100 days, 
~ 30 - 40 percent of the events observed may indeed be related 
to BHs. 

Up to date, thousands of events have been observed toward 
the Galactic bulge by the various surveys (lMoniezil2010h . Thus 
our results suggest that at least several hundreds of events re- 
lated to ICoRs could be present in the catalogs. These are likely 
to be hidden among long-duration events. An exces s of long- 
duratio n events has indeed been reported by iPopowski et al.l 
(I2005h . These event are likely to be generated by massive ob- 
jects, like NSs and BHs. 

The basic problem is therefore to discuss a procedure to 
distinguish isolated compact remnants from normal stars that 
are responsible for the microlensing events. While this discus- 
sion is out of the scope of the present paper, it is worth to 
add some considerations. The only source of steady luminos- 
ity for ICoRs suggested thus far is accretion from the ISM 
(e.glOstriker et aljil970l iTreves & Coi plll99lL iBl aes & Madaul 
I1993L lAeol & Kamionkowskill2002l) . However, in the case of 
NSs, the absence of sources in the ac creting phase as clearly 
demo nstrated by the ROSAT surveys (Neuhauser & Triimper, 
1999) indicates that accretion is inhibited most pr obably 
by the rotating m agnetic structure (e.g. Toro pina et al.l |2003, 
iPerna et aUl2003l) . Accretion onto isolated BHs is favored with 
respect to the NSs, because they are more massive and do 
not have a magnetic field that can hamper the accretion flow. 
However, they are far less numerous than NSs. 

These considerations should be revisited, because we now 
consider the remnants in the bulge. The magnetic field of NSs 
and the properties of the ISM are different from what was dis- 
cussed so far (e.g. IZane et alJl l996). One should first calculate a 
realis tic X-ray luminosity (see e.g. on this line Bold in & Popovl 
2010 and references therein), corrected for absorption, and then 
compare with present and future X-ray missions. The advan- 
tage with respect to a blind search of ICoRs is that microlensing 
events give a precise location of the object. 

We note that, if our calculations are correct, microlensing 
could be the only way to probe the velocity distribution of iso- 
lated BHs, which has not yet been constrained. 



4. Discussion 

Stimulated by recent results on the distribution of ONS in the 
Milky Way, we estimated the contribution of isolated compact 
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Fig. 3. Optical depth profiles of normal stars (left) and ICoRs (right) as a function of the Galactic coordinates. Contour profiles are 
drawn at (0.5, 1.0, 2.0, 3.0) x lO" 6 and (0.5, 1.0, 3.0, 6.0, 8.0) x 10~ 8 . Darker gray levels indicate larger optical depth. 
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Table 3. Rates and average time scales toward the Galactic center and the Baade's window. 
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Table 4. Fractional contribution of NSs and BHs to the rate of events. The l.o.s selected are the same as in Table . 
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Fig. 4. Distribution of time scales for two different l.o.s., (/, b) = (0°, 0°) (left) and (/, b) = (1°, -3°.9) (right). Each panel shows the 
total contribution (solid) and those of normal stars (dotted), NSs (dashed) and BHs (dot-dashed), respectively. 
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Fig. 5. Fractional contribution of ICoRs to the microlensing rate for NSs (left panels) and BHs (right panels). The dashed line 
represents the expected contribution when kinematic effects are not taken into account. 
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